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Summary. I discuss the use of Type la supernovae (SNe la) for cosmological dis- 
tance determinations. Low-redshift SNe la (z ^ 0.1) demonstrate that the Hubble 
expansion is linear, that Ho — 72 ± 8 km s _1 Mpc -1 , and that the properties of 
dust in other galaxies are similar to those of dust in the Milky Way. The light 
curves of high-redshift (z = 0.3-1) SNe la are stretched in a manner consistent 
with the expansion of space; similarly, their spectra exhibit slower temporal evolu- 
tion (by a factor of 1 + z) than those of nearby SNe la. The measured luminosity 
distances of SNe la as a function of redshift have shown that the expansion of the 
Universe is currently accelerating, probably due to the presence of repulsive dark 
energy such as Einstein's cosmological constant (A). From about 200 SNe la, we 
find that H t = 0.96 ± 0.04, and Q A - \AQ M = 0.35 ± 0.14. Combining our data 
with the results of large-scale structure surveys, we find a best fit for Qm and Qa 
of 0.28 and 0.72, respectively — essentially identical to the recent WMAP results 
(and having comparable precision). The sum of the densities, ~ 1.0, agrees with 
extensive measurements of the cosmic microwave background radiation, including 
WMAP, and coincides with the value predicted by most inflationary models for the 
early Universe: the Universe is flat on large scales. A number of possible system- 
atic effects (dust, supernova evolution) thus far do not seem to eliminate the need 
for Qa > 0. However, during the past few years some very peculiar low-redshift 
SNe la have been discovered, and we must be mindful of possible systematic effects 
if such objects are more abundant at high redshifts. Recently, analyses of SNe la 
at z — 1.0-1.7 provide further support for current acceleration, and give tentative 
evidence for an early epoch of deceleration. The dynamical age of the Universe is 
estimated to be 13.1±1.5 Gyr, consistent with the ages of globular star clusters and 
with the WMAP result of 13.7 ± 0.2 Gyr. Current projects include the search for 
additional SNe la at z > 1 to confirm the early deceleration, and the measurement 
of a few hundred SNe la at z — 0.2-0.8 to more accurately determine the equation 
of state of the dark energy, w — P/(pc 2 ), whose value is now constrained by SNe la 
to be in the range —1.48 ,$ w £ — 0.72 at 95% confidence. 



1 Introduction 

Supernovae (SNe) come in two main varieties (see Filippenko 1997b for a 
review) . Those whose optical spectra exhibit hydrogen are classified as Type 
II, while hydrogen-deficient SNe are designated Type I. SNe I are further 
subdivided according to the appearance of the early-time spectrum: SNe la 
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are characterized by strong absorption near 6150 A ( now attributed to Si II) , 
SNe lb lack this feature but instead show prominent He I lines, and SNe Ic 
have neither the Si II nor the He I lines. SNe la are believed to result from 
the thermonuclear disruption of carbon-oxygen white dwarfs, while SNe II 
come from core collapse in massive supergiant stars. The latter mechanism 
probably produces most SNe Ib/Ic as well, but the progenitor stars previously 
lost their outer layers of hydrogen or even helium. 

It has long been recognized that SNe la may be very useful distance 
indicators for a number of reasons; see Branch & Tammann (1992), Branch 
(1998), and references therein. (1) They are exceedingly luminous, with peak 
M B averaging -19.0 mag if H = 72 km s" 1 Mpc" 1 . (2) "Normal" SNe la 
have small dispersion among their peak absolute magnitudes (a ^ 0.3 mag). 
(3) Our understanding of the progenitors and explosion mechanism of SNe la 
is on a reasonably firm physical basis. (4) Little cosmic evolution is expected 
in the peak luminosities of SNe la, and it can be modeled. This makes SNe la 
superior to galaxies as distance indicators. (5) One can perform local tests of 
various possible complications and evolutionary effects by comparing nearby 
SNe la in different environments. 

Research on SNe la in the 1990s has demonstrated their enormous poten- 
tial as cosmological distance indicators. Although there are subtle effects that 
must indeed be taken into account, it appears that SNe la provide among the 
most accurate values of H , qo (the deceleration parameter), Qm (the matter 
density), and Qa [the cosmological constant, Ac 2 / (3Hq)}. 

There have been two major teams involved in the systematic investi- 
gation of high-rcdshift SNe la for cosmological purposes. The "Supernova 
Cosmology Project" (SCP) is led by Saul Perlmutter of the Lawrence Berke- 
ley Laboratory, while the "High-Z Supernova Search Team" (HZT) is led by 
Brian Schmidt of the Mt. Stromlo and Siding Springs Observatories. I have 
been privileged to work with both teams (see Filippenko 2001 for a personal 
account), but my primary allegiance is now with the HZT. 

2 Homogeneity and Heterogeneity 

Until the mid-1990s, the traditional way in which SNe la were used for cosmo- 
logical distance determinations was to assume that they are perfect "standard 
candles" and to compare their observed peak brightness with those of SNe la 
in galaxies whose distances had been independently determined (e.g., with 
Cepheid variables). The rationale was that SNe la exhibit relatively little 
scatter in their peak blue luminosity (<tb ~ 0.4-0.5 mag; Branch & Miller 
1993), and even less if "peculiar" or highly reddened objects were eliminated 
from consideration by using a color cut. Moreover, the optical spectra of 
SNe la are usually rather homogeneous, if care is taken to compare objects 
at similar times relative to maximum brightness (Riess et al. 1997, and ref- 
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erences therein). Over 80% of all SNe la discovered through the early 1990s 
were "normal" (Branch, Fisher, & Nugent 1993). 

From a Hubble diagram constructed with unreddened, moderately distant 
SNe la (z £ 0.1) for which peculiar motions are small and relative distances 
(given by ratios of redshifts) are accurate, Vaughan et al. (1995) find that 

(M B (max)) = (-19.74 ± 0.06) + 5 log (H o /50) mag. (1) 

In a series of papers, Sandage et al. (1996) and Saha et al. (1997) combine 
similar relations with Hubble Space Telescope (HST) Cepheid distances to 
the host galaxies of seven SNe la to derive H — 57 ± 4 km s _1 Mpc -1 . 

Over the past two decades it has become clear, however, that SNe la do 
not constitute a perfectly homogeneous subclass (e.g., Filippcnko 1997a, b). 
In retrospect this should have been obvious: the Hubble diagram for SNe la 
exhibits scatter larger than the photometric errors, the dispersion actually 
rises when reddening corrections are applied (under the assumption that all 
SNe la have uniform, very blue intrinsic colors at maximum; van den Bergh & 
Pazder 1992; Sandage & Tammann 1993), and there are some significant out- 
liers whose anomalous magnitudes cannot be explained by extinction alone. 

Spectroscopic and photometric peculiarities have been noted with increas- 
ing frequency in well-observed SNe la. A striking case is SN 1991T; its pre- 
maximum spectrum did not exhibit Si II or Ca II absorption lines, yet two 
months past maximum brightness the spectrum was nearly indistinguishable 
from that of a classical SN la (Filippenko et al. 1992b; Phillips et al. 1993). 
The light curves of SN 1991T were slightly broader than the SN la tem- 
plate curves, and the object was probably somewhat more luminous than 
average at maximum. Another well-observed, peculiar SNe la is SN 1991bg 
(Filippenko et al. 1992a; Leibundgut et al. 1993; Turatto et al. 1996). At 
maximum brightness it was subluminous by 1.6 mag in V and 2.5 mag in B, 
its colors were intrinsically red, and its spectrum was peculiar (with a deep 
absorption trough due to Ti II). Moreover, the decline from maximum was 
very steep, the /-band light curve did not exhibit a secondary maximum like 
normal SNe la, and the velocity of the ejecta was unusually low. The photo- 
metric heterogeneity among SNe la is well demonstrated by Suntzeff (1996) 
with objects having excellent BVRI light curves. 

3 Cosmological Uses: Low Redshifts 

Although SNe la can no longer be considered perfect "standard candles," 
they are still exceptionally useful for cosmological distance determinations. 
Excluding those of low luminosity (which are hard to find, especially at large 
distances), most of the nearby SNe la that had been discovered through the 
early 1990s were nearly standard (Branch et al. 1993; but see Li et al. 2001b 
for more recent evidence of a higher intrinsic peculiarity rate). Also, after 
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many tenuous suggestions (e.g., Pskovskii 1977, 1984; Branch 1981), Phillips 
(1993) found convincing evidence for a correlation between light curve shape 
and the luminosity at maximum brightness by quantifying the photomet- 
ric differences among a set of nine well-observed SNe la, using a parameter 
[Ami^(B)] that measures the total drop (in B magnitudes) from _B-band 
maximum to t = 15 days later. In all cases the host galaxies of his SNe la 
have accurate relative distances from surface brightness fluctuations or from 
the Tully-Fishcr relation. The intrinsically bright SNe la clearly decline more 
slowly than dim ones, but the correlation is stronger in B than in V or /. 

Using SNe la discovered during the Calan/Tololo survey (z ^ 0.1), Hamuy 
et al. (1995, 1996b) refine the Phillips (1993) correlation between peak lumi- 
nosity and Am\${B). Apparently the slope is steep only at low luminosities; 
thus, objects such as SN 1991bg skew the slope of the best-fitting single 
straight line. Hamuy et al. reduce the scatter in the Hubble diagram of nor- 
mal, unreddened SNe la to only 0.17 mag in B and 0.14 mag in V; see also 
Tripp (1997). Yet another parameterization is the "stretch" method of Perl- 
mutter et al. (1997) and Goldhaber et al. (2001): the i?-band light curves of 
SNe la appear nearly identical when expanded or contracted temporally by 
a factor (1 + s), where the value of s varies among objects. In a similar but 
distinct effort, Riess, Press, & Kirshner (1995) show that the luminosity of 
SNe la correlates with the detailed shape of the overall light curve. 

By using light curve shapes measured through several different filters, 
Riess, Press, & Kirshner (1996a) extend their analysis and objectively elimi- 
nate the effects of interstellar extinction: a SN la that has an unusually red 
B — V color at maximum brightness is assumed to be intrinsically sublumi- 
nous if its light curves rise and decline quickly, or of normal luminosity but 
significantly reddened if its light curves rise and decline more slowly. With a 
set of 20 SNe la consisting of the Calan/Tololo sample and their own objects, 
they show that the dispersion decreases from 0.52 mag to 0.12 mag after ap- 
plication of this "multi-color light curve shape" (MLCS) method. The results 
from an expanded set of nearly 50 SNe la indicate that the dispersion de- 
creases from 0.44 mag to 0.15 mag (Figure 1). The resulting Hubble constant 
is 65 ± 2 (statistical) ±7 (systematic) km s _1 Mpc -1 , with an additional 
systematic and zeropoint uncertainty of ±5 km s -1 Mpc -1 . (Re-calibrations 
of the Cepheid distance scale, and other recent refinements, lead to a best 
estimate of Ho = 72 ± 8 km s -1 Mpc -1 , where the error bar includes both 
statistical and systematic uncertainties; Parodi et al. 2000; Freeman et al. 
2001.) Riess et al. (1996a) also show that the Hubble flow is remarkably lin- 
ear; indeed, SNe la now constitute the best evidence for linearity. Finally, 
they argue that the dust affecting SNe la is not of circumstellar origin, and 
show quantitatively that the extinction curve in external galaxies typically 
does not differ from that in the Milky Way (cf. Branch & Tammann 1992, 
but see Tripp 1998). 
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Figure 1: Hubble diagrams for SNe la (A. G. Ricss 2001, private commu- 
nication) with velocities (km s _1 ) in the COBE rest frame on the Cepheid 
distance scale. The ordinate shows distance modulus, m — M (mag). Top: 
The objects are assumed to be standard candles and there is no correction 
for extinction; the result is a = 0.42 mag and Ho — 58 ± 8 km s _1 Mpc -1 . 
Bottom: The same objects, after correction for reddening and intrinsic differ- 
ences in luminosity. The result is a = 0.15 mag and Ho = 65 ± 2 (statistical) 
km s" 1 Mpc -1 , subject to changes in the zeropoint of the Cepheid distance 
scale. (Indeed, the latest results with SNe la favor H = 72 km s -1 Mpc -1 .) 

The advantage of systematically correcting the luminosities of SNe la 
at high redshifts rather than trying to isolate "normal" ones seems clear in 
view of evidence that the luminosity of SNe la may be a function of stellar 
population. If the most luminous SNe la occur in young stellar populations 
(e.g., Hamuy et al. 1996a, 2000; Branch, Baron, & Romanishin 1996; Ivanov, 
Hamuy, & Pinto 2000), then we might expect the mean peak luminosity of 
high-z SNe la to differ from that of a local sample. Alternatively, the use 
of Cepheids (Population I objects) to calibrate local SNe la can lead to a 
zeropoint that is too luminous. On the other hand, as long as the physics 
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of SNe la is essentially the same in young stellar populations locally and at 
high redshift, we should be able to adopt the luminosity correction methods 
(photometric and spectroscopic) found from detailed studies of low- z SNe la. 

Large numbers of nearby SNe la are now being found by my team's Lick 
Observatory Supernova Search (LOSS) conducted with the 0.76-m Katzman 
Automatic Imaging Telescope (KAIT; Li et al. 2000; Filippenko et al. 2001; 
Filippenko 2003; see http://astro.berkeley.edu/~bait/kait.html). CCD im- 
ages are taken of ~ 1000 galaxies per night and compared with KAIT "tem- 
plate images" obtained earlier; the templates are automatically subtracted 
from the new images and analyzed with computer software. The system re- 
observes the best candidates the same night, to eliminate star-like cosmic 
rays, asteroids, and other sources of false alarms. The next day, undergrad- 
uate students at UC Berkeley examine all candidates, including weak ones, 
and they glance at all subtracted images to locate SNe that might be near 
bright, poorly subtracted stars or galactic nuclei. LOSS discovered 20 SNe 
(of all types) in 1998, 40 in 1999, 38 in 2000, 68 in 2001, and 82 in 2002, mak- 
ing it by far the world's leading search for nearby SNe. The most important 
objects were photometrically monitored through BVRI (and sometimes U) 
filters (e.g., Li et al. 2001a, 2003; Modjaz et al. 2001; Leonard et al. 2002a,b), 
and unfiltered follow-up observations (e.g., Matheson et al. 2001) were made 
of most of them during the course of the SN search. This growing sample of 
well-observed SNe la should allow us to more precisely calibrate the MLCS 
method, as well as to look for correlations between the observed properties 
of the SNe and their environment (Hubble type of host galaxy, metallicity, 
stellar population, etc.). 

4 Cosmological Uses: High Redshifts 

These same techniques can be applied to construct a Hubble diagram with 
high-redshift SNe la, from which the value of qo = (J?m/2) — ft a can be 
determined. With enough objects spanning a range of redshifts, we can mea- 
sure Qm and ft a independently (e.g., Goobar & Perlmutter 1995). Contours 
of peak apparent i?-band magnitude for SNe la at two redshifts have differ- 
ent slopes in the {2m~Qa plane, and the regions of intersection provide the 
answers we seek. 

4.1 The Search 

Based on the pioneering work of Norgaard-Nielsen et al. (1989), whose goal 
was to find SNe in moderate-redshift clusters of galaxies, the SCP (Perlmutter 
et al. 1995a, 1997) and the HZT (Schmidt et al. 1998) devised a strategy that 
almost guarantees the discovery of many faint, distant SNe la "on demand," 
during a predetermined set of nights. This "batch" approach to studying 
distant SNe allows follow-up spectroscopy and photometry to be scheduled in 
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advance, resulting in a systematic study not possible with random discoveries. 
Most of the searched fields are equatorial, permitting follow-up from both 
hemispheres. The SCP was the first group to convincingly demonstrate the 
ability to find SNe in batches. 

Our approach is simple in principle; see Schmidt et al. (1998) for details, 
and for a description of our first high-redshift SN la (SN 1995K). Pairs of first- 
epoch images are obtained with wide-field cameras on large telescopes (e.g., 
the Big Throughput Camera on the CTIO 4-m Blanco telescope) during the 
nights around new moon, followed by second-epoch images 3-4 weeks later. 
(Pairs of images permit removal of cosmic rays, asteroids, and distant Kuiper- 
belt objects.) These are compared immediately using well-tested software, 
and new SN candidates are identified in the second-epoch images (Figure 
2). Spectra are obtained as soon as possible after discovery to verify that 
the objects are SNe la and determine their rcdshifts. Each team has already 
found over nearly 200 SNe in concentrated batches, as reported in numerous 
IAU Circulars (e.g., Perlmutter et al. 1995b, 11 SNe with 0.16 < z < 0.65; 
Suntzcff et al. 1996, 17 SNe with 0.09 < z < 0.84). The observed SN la rate 
at z w 0.5 is consistent with the low-z SN la rate together with plausible 
star- formation histories (Pain et al. 2002; Tonry et al. 2003), but the error 
bars on the high-z rate are still quite large. 




Figure 2: Discovery image of SN 1997cj (28 April 1997), along with the tem- 
plate image and an HST image obtained subsequently. The net (subtracted) 
image is also shown. 
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Intensive photometry of the SNe la commences within a few days after 
procurement of the second-epoch images; it is continued throughout the en- 
suing and subsequent dark runs. In a few cases HST images are obtained. As 
expected, most of the discoveries are on the rise or near maximum bright- 
ness. When possible, the SNe are observed in filters that closely match the 
redshifted B and V bands; this way, the K-corrections become only a second- 
order effect (Kim, Goobar, & Perlmutter 1996; Nugent, Kim, & Perlmutter 
2002). We try to obtain excellent multi-color light curves, so that reddening 
and luminosity corrections can be applied (Riess et al. 1996a; Hamuy et al. 
1996a,b). 

Although SNe in the magnitude range 22-22.5 can sometimes be spec- 
troscopically confirmed with 4-m class telescopes, the signal-to-noise ratios 
are low, even after several hours of integration. Certainly Keck, Gemini, the 
VLT, or Magellan are required for the fainter objects (22.5-24.5 mag). With 
the largest telescopes, not only can we rapidly confirm a substantial num- 
ber of candidate SNe, but we can search for peculiarities in the spectra that 
might indicate evolution of SNe la with redshift. Moreover, high-quality spec- 
tra allow us to measure the age of a SN: we have developed a method for 
automatically comparing the spectrum of a SN la with a library of spectra 
from many different epochs in the development of SNe la (Riess et al. 1997). 
Our technique also has great practical utility at the telescope: we can de- 
termine the age of a SN "on the fly," within half an hour after obtaining its 
spectrum. This allows us to decide rapidly which SNe are best for subsequent 
photometric follow-up, and we immediately alert our collaborators elsewhere. 

4.2 Results 

First, we note that the light curves of high-redshift SNe la are broader than 
those of nearby SNe la; the initial indications (Lcibundgut et al. 1996; Gold- 
haber et al. 1997), based on small numbers of SNe la, are amply confirmed 
with the larger samples (Goldhaber et al. 2001). Quantitatively, the amount 
by which the light curves are "stretched" is consistent with a factor of 1 + z, 
as expected if redshifts are produced by the expansion of space rather than by 
"tired light" and other non-expansion hypotheses for the redshifts of objects 
at cosmological distances. [For non-standard cosmological interpretations of 
the SN la data, see Narlikar & Arp (1997) and Hoyle, Burbidge, & Narlikar 
(2000).] We also demonstrate this spectroscopically at the 2cr confidence level 
for a single object: the spectrum of SN 1996bj (z = 0.57) evolved more slowly 
than those of nearby SNe la, by a factor consistent with 1 + z (Riess et al. 
1997). More recently, we have used observations of SN 1997ex (z = 0.36) 
at three epochs to conclusively verify the effects of time dilation: temporal 
changes in the spectra are slower than those of nearby SNe la by roughly the 
expected factor of 1.36. Although one might be able to argue that something 
other than universal expansion could be the cause of the apparent stretching 
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of SN la light curves at high redshifts, it is much more difficult to attribute 
apparently slower evolution of spectral details to an unknown effect. 

The formal value of Qm derived from SNe la has changed with time. 
The SCP published the first result (Perlmutter et al. f995a), based on a 
single object, SN 1992bi at z = 0.458: Q M = 0.2 ± 0.6 ± 1.1 (assuming that 
Qa = 0). The SCP's analysis of their first seven objects (Perlmutter et al. 
1997) suggested a much larger value of Qm = 0.88 ± 0.6 (if Qa = 0) or 
Qm = 0.94 ±0.3 (if J?totai = !)• Such a high-density universe seemed at odds 
with other, independent measurements of Qm- However, with the subsequent 
inclusion of just one more object, SN 1997ap at z — 0.83 (the highest known 
for a SN la at the time; Perlmutter et al. 1998), their estimates were revised 
back down to Qm — 0.2 ±0.4 if Qa = 0, and Qm = 0.6 ±0.2 if J?totai = 1; the 
apparent brightness of SN 1997ap had been precisely measured with HST, so 
it substantially affected the best fits. 

Meanwhile, the HZT published (Garnavich et al. 1998a) an analysis of 
four objects (three of them observed with HST), including SN 1997ck at 
z = 0.97, at that time a redshift record, although they cannot be absolutely 
certain that the object was a SN la because the spectrum is too poor. From 
these data, the HZT derived that Qm = —0.1 ±0.5 (assuming Qa = 0) 
and Qm = 0.35 ± 0.3 (assuming J?totai = 1), inconsistent with the high Qm 
initially found by Perlmutter et al. (1997) but consistent with the revised 
estimate in Perlmutter et al. (1998). An independent analysis of 10 SNe la 
using the "snapshot" distance method (with which conclusions are drawn 
from sparsely observed SNe la) gave quantitatively similar conclusions (Ricss 
et al. 1998a). However, none of these early data sets carried the statistical 
discriminating power to detect cosmic acceleration. 

The SCP's next results were announced at the 1998 January AAS meet- 
ing in Washington, DC. A press conference was scheduled, with the stated 
purpose of presenting and discussing the then-current evidence for a \ow-Qm 
universe as published by Perlmutter et al. (1998; SCP) and Garnavich et al. 
(1998a; HZT). When showing the SCP's Hubble diagram for SNe la, however, 
Perlmutter also pointed out tentative evidence for acceleration*. He said that 
the conclusion was uncertain, and that the data were equally consistent with 
no acceleration; the systematic errors had not yet been adequately assessed. 
Essentially the same conclusion was given by the SCP in their talks at a 
conference on dark matter, near Los Angeles, in February 1998 (Goldhaber 
& Perlmutter 1998). 

Although it chose not to reveal them at the same 1998 January AAS 
meeting, the HZT already had similar, tentative evidence for acceleration in 
their own SN la data set. The HZT continued to perform numerous checks of 
their data analysis and interpretation, including fairly thorough consideration 
of various possible systematic effects. Unable to find any significant problems, 
even with the possible systematic effects, the HZT reported detection of a 
nonzero value for Qa (based on 16 high-z SNe la) at the Los Angeles dark 
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matter conference in February 1998 (Filippenko & Riess 1998), and soon 
thereafter submitted a formal paper that was published in September 1998 
(Riess ct al. 1998b). Their original Hubble diagram for the 10 best-observed 
high-z SNe la is given in Figure 3. With the MLCS method applied to the 
full set of 16 SNe la, the HZT's formal results were Q M = 0.24 ± 0.10 if 
Aotai = 1, or J?m = —0.35 ± 0.18 (unphysical) if Qa = 0. If one demanded 
that Qm — 0.2, then the best value for Qa was 0.66±0.21. These conclusions 
did not change significantly when only the 10 best-observed SNe la were used 
(Figure 3; Q M = 0.28 ± 0.10 if Q totaX = 1). 




z 



Figure 3 (left): The upper panel shows the Hubble diagram for the low-z and 
high-z HZT SN la sample with MLCS distances; see Riess et al. (1998b). 
Overplotted are three world models: "low" and "high" Qm with Qa = 0, and 
the best fit for a flat universe {Qm = 0.28, Qa = 0.72). The bottom panel 
shows the difference between data and models from the Qm = 0.20, Qa = 
prediction. Only the 10 best-observed high-z SNe la are shown. The average 
difference between the data and the Qm = 0.20, Qa — prediction is ~ 0.25 
mag. 

Figure 4 (right): The HZT's combined constraints from SNe la (left) and 
the position of the first acoustic peak of the cosmic microwave background 
(CMB) angular power spectrum, based on data available in mid-1998; see 
Garnavich et al. (1998b). The contours mark the 68.3%, 95.4%, and 99.7% en- 
closed probability regions. Solid curves correspond to results from the MLCS 
method; dotted ones are from the Ami$(B) method; all 16 SNe la in Riess 
et al. (1998b) were used. 
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Another important constraint on the cosmological parameters could be 
obtained from measurements of the angular scale of the first acoustic peak 
of the CMB (e.g., Zaldarriaga, Spergel, & Seljak 1997; Eisenstein, Hu, & 
Tegmark 1998); the SN la and CMB techniques provide nearly complemen- 
tary information. A stunning result was already available by mid-1998 from 
existing measurements (e.g., Hancock et al. 1998; Lineweaver & Barbosa 
1998): the HZT's analysis of the SN la data in Riess et al. (1998b) demon- 
strated that fl M + H A = 0.94 ± 0.26 (Figure 4), when the SN and CMB 
constraints were combined (Garnavich et al. 1998b; see also Lineweaver 1998, 
Efstathiou et al. 1999, and others). 

Somewhat later (June 1999), the SCP published almost identical results, 
implying an accelerating expansion of the Universe, based on an essentially 
independent set of 42 high-z SNc la (Perlmutter et al. 1999). Their data, 
together with those of the HZT, are shown in Figure 5, and the corresponding 
confidence contours in the Qa vs. J?m plane are given in Figure 6. This 
incredible agreement suggested that neither group had made a large, simple 
blunder; if the result was wrong, the reason must be subtle. Had there been 
only one team working in this area, it is likely that far fewer astronomers and 
physicists throughout the world would have taken the result seriously. 

Moreover, already in 1998-1999 there was tentative evidence that the 
"dark energy" driving the accelerated expansion was indeed consistent with 
the cosmological constant, A. If A dominates, then the equation of state of 
the dark energy should have an index w = — 1, where the pressure (P) and 
density (p) are related according to w = P/(pc 2 ). Garnavich et al. (1998b) 
and Perlmutter et al. (1999) were able to set an interesting limit, w ;$ — 0.60 
at the 95% confidence level. However, more high-quality data at z « 0.5 
are needed to narrow the allowed range, in order to test other proposed 
candidates for dark energy such as various forms of "quintessence" (e.g., 
Caldwell, Dave, & Steinhardt 1998). 

Although the CMB results appeared reasonably persuasive in 1998-1999, 
one could argue that fluctuations on different scales had been measured with 
different instruments, and that suble systematic effects might lead to erro- 
neous conclusions. These fears were dispelled only 1-2 years later, when the 
more accurate and precise results of the BOOMERANG collaboration were 
announced (do Bernardis et al. 2000, 2002). Shortly thereafter the MAXIMA 
collaboration distributed their very similar findings (Hanany et al. 2000; Balbi 
et al. 2000; Netterfield et al. 2002; see also the TOCO, DASI, and many other 
measurements) . Figure 6 illustrates that the CMB measurements tightly con- 
strain J?totai t° be close to unity; we appear to live in a flat universe, in agree- 
ment with most inflationary models for the early Universe! Combined with 
the SN la results, the evidence for nonzero Qa was fairly strong. Making the 
argument even more compelling was the fact that various studies of clusters 
of galaxies (see summary by Bahcall et al. 1999) showed that Qm ~ 0.3, 
consistent with the results in Figures 4 and 6. Thus, a "concordance cosmol- 
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ogy" had emerged: Qm ~ 0.3, « 0.7 — consistent with what had been 
suspected some years earlier by Ostriker & Steinhardt (1995; see also Carroll, 
Press, & Turner 1992). 
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Figure 5 (left): As in Figure 3, but this time including both the HZT (Riess 
et al. 1998b) and SCP (Perlmutter et al. 1999) samples of low-redshift and 
high-rcdshift SNe la. Overplotted are three world models: Qm = 0.3 and 1.0 
with = 0, and a flat universe (/2totaJ = 1*0) with Qa = 0.7. The bottom 
panel shows the difference between data and models from the Dm = 0.3, 
Qa — prediction. 

Figure 6 (right): The combined constraints from SNe la (see Figure 5) and 
the position of the first acoustic peak of the CMB angular power spectrum, 
based on BOOMERANG and MAXIMA data. The contours mark the 68.3%, 
95.4%, and 99.7% enclosed probability regions determined from the SNe la. 
According to the CMB, f? t otal ~ 1-0. 

Yet another piece of evidence for a nonzero value of A was provided 
by the Two-Degree Field Galaxy Redshift Survey (2dFGRS; Peacock et al. 
2001; Percival et al. 2001; Efstathiou et al. 2002). Combined with the CMB 
maps, their results are inconsistent with a universe dominated by gravitating 
dark matter. Again, the implication is that about 70% of the mass-energy 
density of the Universe consists of some sort of dark energy whose gravita- 
tional effect is repulsive. Very recently, results from the Wilkinson Microwave 
Anisotropy Prove (WMAP) appeared; together with the 2dFGRS constraints, 
they confirm and refine the concordance cosmology (Dm — 0.27, Qa — 0.73, 
j? baryon = 0.044, H a = 71 ± 4 km s" 1 Mpc" 1 ; Spergel et al. 2003). 

The dynamical age of the Universe can be calculated from the cosmolog- 
ical parameters. In an empty Universe with no cosmological constant, the 
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dynamical age is simply the "Hubble time" to (i- e -7 the inverse of the Hubble 
constant); there is no deceleration. In the late-1990s, SNe la gave H — 65 ±7 
km s _1 Mpc -1 , and a Hubble time of 15.1 ± 1.6 Gyr. For a more complex 
cosmology, integrating the velocity of the expansion from the current epoch 
(z = 0) to the beginning (z — oo) yields an expression for the dynamical age. 
As shown in detail by Riess et al. (1998b), by mid-1998 the HZT had obtained 
a value of 14.2+ ; 8 Gyr (with Hq — 65) using the likely range for {Qm, Qa) 
that they measured. (The precision was so high because their experiment was 
sensitive to roughly the difference between Qm and Qa, and the dynamical 
age also varies in approximately this way.) Including the systematic uncer- 
tainty of the Ccpheid distance scale, which may be up to 10%, a reasonable 
estimate of the dynamical age was 14.2 ±1.7 Gyr (Riess et al. 1998b). Again, 
the SCP's result was very similar (Perlmutter et al. 1999), since it was based 
on nearly the same derived values for the cosmological parameters. The most 
recent results, reported by Tonry et al. (2003) and adopting H n = 72 ± 8 km 
s -1 Mpc -1 , give a dynamical age of 13.1 ± 1.5 Gyr for the Universe — again, 
in agreement with the WMAP result of 13.7 ± 0.2 Gyr. 

This expansion age is also consistent with ages determined from various 
other techniques such as the cooling of white dwarfs (Galactic disk > 9.5 
Gyr; Oswalt et al. 1996), radioactive dating of stars via the thorium and 
europium abundances (15.2 ± 3.7 Gyr; Cowan et al. 1997), and studies of 
globular clusters (10-15 Gyr, depending on whether Hipparcos parallaxes of 
Ccpheids are adopted; Gratton et al. 1997; Chaboyer et al. 1998). The ages of 
the oldest stars no longer seem to exceed the expansion age of the Universe; 
the long-standing "age crisis" has evidently been resolved. 

5 Discussion 

Although the convergence of different methods on the same answer is reassur- 
ing, and suggests that the concordance cosmology is correct, it is important 
to vigorously test each method to make sure it is not leading us astray. More- 
over, only through such detailed studies will the accuracy and precision of 
the methods improve, allowing us to eventually set better constraints on the 
equation of state parameter, w. Here I discuss the systematic effects that 
could adversely affect the SN la results. 

High-rcdshift SNe la are observed to be dimmer than expected in an 
empty Universe (i.e., Qm = 0) with no cosmological constant. At z w 0.5, 
where the SN la observations have their greatest leverage on A, the difference 
in apparent magnitude between an Qm = 0.3 {Qa = 0) universe and a flat 
universe with Qa = 0.7 is only about 0.25 mag. Thus, we need to find out if 
chemical abundances, stellar populations, selection bias, gravitational lensing, 
or grey dust can have an effect this large. Although both the HZT and SCP 
had considered many of these potential systematic effects in their original 
discovery papers (Riess et al. 1998b; Perlmutter et al. 1999), and had shown 
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with reasonable confidence that obvious ones were not greatly affecting their 
conclusions, if was of course possible that they were wrong, and that the data 
were being misinterpreted. 

5.1 Evolution 

Perhaps the most obvious possible culprit is evolution of SNe la over cosmic 
time, due to changes in metallicity, progenitor mass, or some other factor. 
If the peak luminosity of SNe la were lower at high redshift, then the case 
for Qa > would weaken. Conversely, if the distant explosions are more 
powerful, then the case for acceleration strengthens. Theorists are not yet sure 
what the sign of the effect will be, if it is present at all; different assumptions 
lead to different conclusions (Hoflich et al. 1998; Umeda et al. 1999; Nomoto 
et al. 2000; Yungelson & Livio 2000). 

Of course, it is extremely difficult, if not effectively impossible, to obtain 
an accurate, independent measure of the peak luminosity of high-z SNe la, 
and hence to directly test for luminosity evolution. However, we can more 
easily determine whether other observable properties of low-z and high-z 
SNe la differ. If they are all the same, it is more probable that the peak 
luminosity is constant as well — but if they differ, then the peak luminosity 
might also be affected (e.g., Hoflich et al. 1998). Drell, Loredo, & Wasserman 
(2000), for example, argue that there are reasons to suspect evolution, because 
the average properties of existing samples of high-z and low-z SNe la seem 
to differ (e.g., the high-z SNe la are more uniform). 

The local sample of SNe la displays a weak correlation between light 
curve shape (or peak luminosity) and host galaxy type, in the sense that the 
most luminous SNe la with the broadest light curves only occur in late-type 
galaxies. Both early-type and late-type galaxies provide hosts for dimmer 
SNe la with narrower light curves (Hamuy et al. 1996a). The mean luminosity 
difference for SNe la in late-type and early-type galaxies is ~ 0.3 mag. In 
addition, the SN la rate per unit luminosity is almost twice as high in late- 
type galaxies as in early-type galaxies at the present epoch (Cappellaro et 
al. 1997). These results may indicate an evolution of SNe la with progenitor 
age. Possibly relevant physical parameters are the mass, metallicity, and C/O 
ratio of the progenitor (Hoflich et al. 1998). 

We expect that the relation between light curve shape and peak luminosity 
that applies to the range of stellar populations and progenitor ages encoun- 
tered in the late-type and early-type hosts in our nearby sample should also 
be applicable to the range we encounter in our distant sample. In fact, the 
range of age for SN la progenitors in the nearby sample is likely to be larger 
than the change in mean progenitor age over the 4-6 Gyr lookback time 
to the high-z sample. Thus, to first order at least, our local sample should 
correct the distances for progenitor or age effects. 
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Figure 7 (left): Spectral comparison (in f x ) of SN 1998ai (z = 0.49; Keck 
spectrum) with low-redshift (z < 0.1) SNe la at a similar age (~ 5 days 
before maximum brightness), from Riess et al. (1998b). The spectra of the 
low-redshift SNe la were resampled and convolved with Gaussian noise to 
match the quality of the spectrum of SN 1998ai. Overall, the agreement in the 
spectra is excellent, tentatively suggesting that distant SNe la are physically 
similar to nearby SNe la. SN 1994B (z = 0.09) differs the most from the 
others, and was included as a "decoy." 

Figure 8 (right): Heavily smoothed spectra of two high-z SNe (SN 1999ff 
at z = 0.455 and SN 1999fv at z = 1.19; quite noisy below -3500 A) are 
presented along with several low-z SN la spectra (SNe 1989B, 1992A, and 
1981B), a SN lb spectrum (SN 1993J), and a SN Ic spectrum (SN 19941); 
see Filippenko (1997) for a discussion of spectra of various types of SNe. The 
date of the spectra relative to -B-band maximum is shown in parentheses 
after each object's name. Specific features seen in SN 1999ff and labeled with 
a letter are discussed by Coil et al. (2000). This comparison shows that the 
two high-z SNe are most likely SNe la. 

We can place empirical constraints on the effect that a change in the 
progenitor age would have on our SN la distances by comparing subsamples 
of low-redshift SNe la believed to arise from old and young progenitors. In the 
nearby sample, the mean difference between the distances for the early-type 
hosts (8 SNe la) and late-type hosts (19 SNe la), at a given redshift, is 0.04 
± 0.07 mag from the MLCS method. This difference is consistent with zero. 
Even if the SN la progenitors evolved from one population at low redshift 
to the other at high redshift, we still would not explain the surplus in mean 
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distance of 0.25 mag over the £Ia — prediction. Moreover, in a major study 
of high-redshift SNe la as a function of galaxy morphology, the SCP found 
no clear differences (except for the amount of scatter; see §5.2) between the 
cosmological results obtained with SNe la in late- type and early-type galaxies 
(Sullivan et al. 2003). 

It is also reassuring that initial comparisons of high-z SN la spectra ap- 
pear remarkably similar to those observed at low redshift. For example, the 
spectral characteristics of SN 1998ai (z = 0.49) appear to be essentially in- 
distinguishable from those of normal low-z SNe la; see Figure 7. In fact, the 
most obviously discrepant spectrum in this figure is the second one from the 
top, that of SN 1994B (z — 0.09); it is intentionally included as a "decoy" 
that illustrates the degree to which even the spectra of nearby, relatively nor- 
mal SNe la can vary. Nevertheless, it is important to note that a dispersion 
in luminosity (perhaps 0.2 mag) exists even among the other, more normal 
SNe la shown in Figure 7; thus, our spectra of SN 1998ai and other high- z 
SNe la are not yet sufficiently good for independent, precise determinations 
of peak luminosity from spectral features (Nugent et al. 1995). Many of them, 
however, are sufficient for ruling out other SN types (Figure 8), or for identi- 
fying gross peculiarities such as those shown by SNe 1991T and 1991bg; see 
Coil et al. (2000). 

We can help verify that the SNe at z ~ 0.5 being used for cosmology do not 
belong to a subluminous population of SNe la by examining restframe /-band 
light curves. Normal, nearby SNe la show a pronounced second maximum in 
the / band about a month after the first maximum and typically about 0.5 
mag fainter (e.g., Ford et al. 1993; Suntzeff 1996). Subluminous SNe la, in 
contrast, do not show this second maximum, but rather follow a linear decline 
or show a muted second maximum (Filippenko et al. 1992a). As discussed by 
Riess et al. (2000), tentative evidence for the second maximum is seen from 
the HZT's existing J-band (restframe /-band) data on SN 1999Q (z = 0.46); 
see Figure 10. Additional tests with spectra and near-infrared light curves 
are currently being conducted. 

Another way of using light curves to test for possible evolution of SNe la 
is to see whether the rise time (from explosion to maximum brightness) is the 
same for high-redshift and low-redshift SNe la; a difference might indicate 
that the peak luminosities are also different (Hoflich et al. 1998). Riess et 
al. (1999c) measured the risetime of nearby SNe la, using data from KAIT, 
the Beijing Astronomical Observatory (BAO) SN search, and a few amateur 
astronomers. Though the exact value of the risetime is a function of peak 
luminosity, for typical low-redshift SNe la it is 20.0 ± 0.2 days. Riess et al. 
(1999b) pointed out that this differs by 5.8<7 from the 'preliminary risetime of 
17.5 ±0.4 days reported in conferences by the SCP (Goldhaber et al. 1998a, b; 
Groom 1998). However, more thorough analyses of the SCP data (Aldcring, 
Knop, & Nugent 2000; Goldhaber et al. 2001) show that the high-redshift 
uncertainty of ±0.4 days that the SCP originally reported was much too small 
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because it did not account for systematic effects. The revised discrepancy with 
the low-redshift risetime is about la or less. Thus, the apparent difference 
in risetimes might be insignificant. Even if the difference is real, however, its 
relevance to the peak luminosity is unclear; the light curves may differ only 
in the first few days after the explosion, and this could be caused by small 
variations in conditions near the outer part of the exploding white dwarf that 
are inconsequential at the peak. 
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Figure 9: The MLCS fit (Riess et al. 1998b; left panel) and the stretch method 
fit (Pcrlmuttcr ct al. 1999; right panel) for SN 2000cx. The MLCS fit is the 
worst we had ever seen through 2000. For the stretch method fit, the solid 
line is the fit to all the data points from t = — 8 to 32 days, the dash-dotted 
line uses only the premaximum datapoints, and the dashed line only the 
postmaximum datapoints. The three fits give very different stretch factors. 
From Li ct al. (2001a). 

Although there are no clear signs that cosmic evolution of SNe la seriously 
compromises our results, it is wise to remain vigilant for possible problems. At 
low redshifts, for example, we already know that some SNe la don't conform 
with the correlation between light curve shape and luminosity. SN 2000cx in 
the SO galaxy NGC 524, for example, has light curves that cannot be fit well 
by any of the fitting techniques currently available (Li et al. 2001a; Filippenko 
2003); see Figure 9. Its late-time color is remarkably blue, inconsistent with 
the homogeneity described by Phillips et al. (1999). The spectral evolution of 
SN 2000cx is peculiar as well: the photosphere appears to have remained hot 
for a long time, and both iron-peak and intermediate-mass elements move 
at very high velocities. An even more peculiar object is SN 2002cx (Li et al. 
2003; Filippenko 2003). It is spectroscopically bizarre, with extremely low 
expansion velocities and almost no evidence for intermediate-mass elements. 
The nebular phase was reached incredibly soon after maximum brightness, 
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despite the low velocity of the ejecta, suggesting that the ejected mass is 
small. SN 2002cx was subluminous by ~ 2 mag at all optical wavelengths 
relative to normal SNe la, despite the early-time spectroscopic resemblance to 
the somewhat over luminous SN 199 IT. The i?-band and /-band light curves 
of SN 2002cx are completely unlike those of normal SNe la. No existing 
theoretical model successfully explains all observed aspects of SN 2002cx. If 
there are more strange beasts like SNe 2000cx and 2002cx at high redshifts 
than at low redshifts, systematic errors may creep into the analysis of high- z 
SNe la. 

5.2 Extinction 

Our SN la distances have the important advantage of including corrections 
for interstellar extinction occurring in the host galaxy and the Milky Way. 
Extinction corrections based on the relation between SN la colors and lumi- 
nosity improve distance precision for a sample of nearby SNe la that includes 
objects with substantial extinction (Riess et al. 1996a); the scatter in the 
Hubble diagram is much reduced. Moreover, the consistency of the measured 
Hubble flow from SNe la with late-type and early-type hosts (see §5.1) shows 
that the extinction corrections applied to dusty SNe la at low redshift do 
not alter the expansion rate from its value measured from SNe la in low-dust 
environments. 

In practice, the high-redshift SNe la generally appear to suffer very little 
extinction; their B — V colors at maximum brightness are normal, suggesting 
little color excess due to reddening. The most detailed available study is that 
of the SCP (Sullivan et al. 2003): they found that the scatter in the Hub- 
ble diagram is minimal for SNe la in early-type host galaxies, but increases 
for SNe la in late-type galaxies. Moreover, on average the SNe in late-type 
galaxies are slightly fainter (by 0.14 ± 0.09 mag) than those in early- type 
galaxies. Finally, at peak brightness the colors of SNe la in late-type galaxies 
are marginally redder than those in early-type galaxies. Sullivan et al. (2003) 
conclude that extinction by dust in the host galaxies of SNe la is one of the 
major sources of scatter in the high-redshift Hubble diagram. By restricting 
their sample to SNe la in early-type host galaxies (presumably with minimal 
extinction), they obtain a very tight Hubble diagram that suggests a nonzero 
value for Qa at the 5<r confidence level, under the assumption that l^totai = 1 • 
In the absence of this assumption, SNe la in early-type hosts still imply that 
Qa > at nearly the 98% confidence level. The results for CLa with SNe la 
in late-type galaxies are quantitatively similar, but statistically less secure 
because of the larger scatter. 

Riess, Press, & Kirshner (1996b) found indications that the Galactic ra- 
tios between selective absorption and color excess are similar for host galaxies 
in the nearby [z < 0.1) Hubble flow. Yet, what if these ratios changed with 
lookback time (e.g., Aguirre 1999a)? Could an evolution in dust-grain size de- 
scending from ancestral interstellar "pebbles" at higher redshifts cause us to 
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underestimate the extinction? Large dust grains would not imprint the red- 
dening signature of typical interstellar extinction upon which our corrections 
necessarily rely. 

However, viewing our SNe through such gray interstellar grains would also 
induce a dispersion in the derived distances. Using the results of Hatano, 
Branch, & Deaton (1998), Riess et al. (1998b) estimate that the expected 
dispersion would be 0.40 mag if the mean gray extinction were 0.25 mag (the 
value required to explain the measured MLCS distances without a cosmo- 
logical constant). This is significantly larger than the 0.21 mag dispersion 
observed in the high-redshift MLCS distances. Furthermore, most of the ob- 
served scatter is already consistent with the estimated statistical errors, leav- 
ing little to be caused by gray extinction. Nevertheless, if we assumed that 
all of the observed scatter were due to gray extinction, the mean shift in the 
SN la distances would be only 0.05 mag. With the existing observations, it 
is difficult to rule out this modest amount of gray interstellar extinction. 
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Figure 10 (left): Restframe /-band (observed J-band) light curve of SN 1999Q 
(z = 0.46, 5 solid points; Riess et al. 2000), and the J-band light curves of 
several nearby SNe la. Subluminous SNe la exhibit a less prominent second 
maximum than do normal SNe la. 

Figure 11 (right): Color excess, Eb-i, for SN 1999Q and different dust models 
(Riess et al. 2000). The data are most consistent with no dust and Qa > 0. 



Gray intergalactic extinction could dim the SNe without either telltale 
reddening or dispersion, if all lines of sight to a given redshift had a similar 
column density of absorbing material. The component of the intergalactic 
medium with such uniform coverage corresponds to the gas clouds producing 
Lyman-a forest absorption at low redshifts. These clouds have individual H I 
column densities less than about 10 15 cm~ 2 (Bahcall et al. 1996). However, 
they display low metallicities, typically less than 10% of solar. Gray extinction 
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would require larger dust grains which would need a larger mass in heavy 
elements than typical interstellar grain size distributions to achieve a given 
extinction. It is possible that large dust grains are blown out of galaxies by 
radiation pressure, and are therefore not associated with Lyman-a clouds 
(Aguirre 1999b). 

But even the dust postulated by Aguirre (1999a, b) and Aguirre & Haiman 
(1999) is not completely gray, having a size of about 0.1 fim. We can test for 
such nearly gray dust by observing high-redshift SNe la over a wide wave- 
length range to measure the color excess it would introduce. If Ay = 0.25 
mag, then E(U-I) and E(B-I) should be 0.12-0.16 mag (Aguirre 1999a,b). 
If, on the other hand, the 0.25 mag faintncss is due to A, then no such red- 
dening should be seen. This effect is measurable using proven techniques; so 
far, with just one SN la (SN 1999Q; Figure 11), our results favor the no-dust 
hypothesis to better than 2a (Riess et al. 2000). More work along these lines 
is in progress. 

5.3 The Smoking Gun 

Suppose, however, that for some reason the dust is very gray, or our color 
measurements are not sufficiently precise to rule out Aguirre's (or other) dust. 
Or, perhaps some other astrophysical systematic effect is fooling us, such as 
possible evolution of the white dwarf progenitors (e.g., Hofiich ct al. 1998; 
Umeda et al. 1999), or gravitational lensing (Wambsganss, Cen, & Ostrikcr 
1998). The most decisive test to distinguish between A and cumulative sys- 
tematic effects is to examine the deviation of the observed peak magnitude 
of SNe la from the magnitude expected in the 1ow-J?m, zero-yl model. If A is 
positive, the deviation should actually begin to decrease at z « 1; we will be 
looking so far back in time that the A effect becomes small compared with 
£2m, and the Universe is decelerating at that epoch. If, on the other hand, a 
systematic bias such as gray dust or evolution of the white dwarf progenitors 
is the culprit, we expect that the deviation of the apparent magnitude will 
continue growing, unless the systematic bias is set up in such an unlikely way 
as to mimic the effects of A (Drcll ct al. 2000). A turnover, or decrease of 
the deviation of apparent magnitude at high redshift, can be considered the 
"smoking gun" of A. 

In a wonderful demonstration of good luck and hard work, Riess et al. 
(2001) report on HST observations of a probable SN la at z w 1.7 (SN 1997ff, 
the most distant SN ever observed) that suggest the expected turnover is 
indeed present, providing a tantalizing glimpse of the epoch of deceleration. 
(See also Bem'tez et al. 2002, which corrects the observed magnitude of SN 
1997ff for gravitational lensing.) SN 1997ff was discovered by Gilliland & 
Phillips (1998) in a repeat HST observation of the Hubble Deep Field-North, 
and serendipitously monitored in the infrared with 7JST/NICMOS. The peak 
apparent SN brightness is consistent with that expected in the decelerating 
phase of the concordance cosmological model, Qm ~ 0.3, Ha ~ 0.7 (Figure 



The Accelerating Universe 21 



12). It is inconsistent with gray dust or simple luminosity evolution, when 
combined with the data for SNe la at z ~ 0.5. On the other hand, it is 
wise to remain cautious: the error bars are large, and it is always possible 
that we are being fooled by this single object. The HZT and SCP currently 
have programs to find and measure more SNe la at such high redshifts. For 
example, SN candidates at very high redshifts (e.g., Giavalisco et al. 2002) 
have been found by "piggybacking" on the Great Observatories Origins Deep 
Survey (GOODS) being conducted with the Advanced Camera for Surveys 
aboard HST; see Figure 13. 
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Figure 12: Hubble diagram for SNe la relative to an empty universe (J? = 0) 
compared with cosmological and astrophysical models (Riess et al. 2001). 
Low-redshift SNe la are from Hamuy et al. (1996a) and Riess et al. (1999a). 
The magnitude of SN 1997ff at z = 1.7 has been corrected for gravitational 
lensing (Bem'tez et al. 2002). The measurements of SN 1997ff are inconsistent 
with astrophysical effects that could mimic previous evidence for an acceler- 
ating universe from SNe la at z s=s 0.5. 
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Figure 13: SN 2002hp, a high-redshift supernova from the GOODS program 
using HST. One can see it brightening and subsequently fading with time. 
The assumed host galaxy is at the top of each frame. 
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Less ambitious programs, concentrating on SNe la at z <L 0.8, have already 
been completed (HZT: Tonry et al. 2003) or are nearing completion (SCP). 
Tonry et al. (2003) measured several SNe la at z w 1, and their deviation 
of apparent magnitude from the \ow-Qm, zero-Tl model is roughly the same 
as that at z « 0.5 (Figure 14), in agreement with expectations based on 
the results of Ricss et al. (2001). Moreover, the new sample of high-redshift 
SNe la presented by Tonry et al., analyzed with methods distinct from (but 
similar to) those used previously, confirm the result of Ricss et al. (1998b) and 
Perlmuttcr et al. (1999) that the expansion of the Universe is accelerating. 




Figure 1^: The new Tonry et al. (2003) data and other data points are shown 
in a residual Hubble diagram: apparent magnitude difference between the 
expected magnitude in an empty universe and the observed magnitude of 
SNe la at each redshift. The highlighted points correspond to median values 
in eight redshift bins. From top to bottom the curves show {Qm,^a) = 
(0.3,0.7), (0.3,0.0), and (1.0,0.0), respectively. 

Figure 15: From Tonry et al. (2003), the probability contours for fl\ versus 
flu are shown at la, 2cr, and 3a [assuming w — —1, where w = P/(pc 2 )]. 
Also shown are the corresponding contours when a prior of fluh = 0.20±0.03 
(where H = lOOh km s" 1 Mpc' 1 ) is adopted from the 2dFGRS (Pcrcival et 
al. 2001). These constraints use the full sample of 172 SNe la with z > 0.01 
and Ay < 0.5 mag. 

By combining all of the available data sets, Tonry et al. (2003) are able 
to use about 200 SNe la, obtaining an incredibly firm detection of Qa > 
(Figure 15). They place the following constraints on cosmological quantities: 
(1) If the equation of state parameter of the dark energy is w = —1, then 
H t a = 0.96 ± 0.04, and Q A - IAQ M = 0.35 ± 0.14. (2) Including the con- 
straint of a flat universe, they find that flu = 0.28±0.05, independent of any 
large-scale structure measurements. (3) Adopting a prior based on the 2dF- 
GRS constraint on flu (Percival et al. 2001) and assuming a flat universe, 
they derive that —1.48 < w < —0.72 at 95% confidence. (4) Adopting the 
2dFGRS results, they find flu = 0.28 and Qa — 0.72, independent of any 



The Accelerating Universe 



23 



assumptions about J?totai- These constraints are similar in precision and in 
value to very recent conclusions reported using WMAP (Spergel et al. 2003), 
also in combination with the 2dFGRS. Complete details on the SN la results 
can be found in Tonry et al. (2003). 

5.4 Measuring the Dark Energy Equation of State 

Every energy component in the Universe can be parameterized by the way its 
density varies as the Universe expands (scale factor a), with p oc aT^ 1+w \ 
and w reflects the component's equation of state, w = P/(pc 2 ), where P 
is the pressure exerted by the component. So for matter, w = 0, while an 
energy component that does not vary with scale factor has w = — 1, as in the 
cosmological constant A. Some really strange energies may have w < — 1: their 
density increases with time (Carroll, Hoffman, & Trodden 2003)! Clearly, a 
good estimate of w becomes the key to differentiating between models. 

The CMB observations imply that the geometry of the universe is close 
to flat, so the energy density of the dark component is simply related to the 
matter density by Q x — 1 — Qm- This allows the luminosity distance as a 
function of rcdshift to be written as 



showing that the dark energy density and equation of state directly influence 
the apparent brightness of standard candles. As demonstrated graphically in 
Figure 16, SNe la observed over a wide range of redshifts can constrain the 
dark energy parameters to a cosmologically interesting accuracy. 

But there are two major problems with using SNe la to measure w. First, 
systematic uncertainties in SN la peak luminosity limit how well P>l{z) can be 
measured. While statistical uncertainty can be arbitrarily reduced by finding 
thousands of SNe la, intrinsic SN properties such as evolution and progen- 
itor mctallicity, and observational limits like photometric calibrations and 
K-corrections, create a systematic floor that cannot be decreased by sheer 
force of numbers. We expect that systematics can be controlled to at best 
3%, with considerable effort. 

Second, SNe at z > 1.0 are very hard to discover and study from the 
ground. As discussed above, both the HZT and the SCP have found a few 
SNe la at z > 1.0, but the numbers and quality of these light curves are 
insufficient for a w measurement. Large numbers of SNe la at z > 1.0 are 
best left to a wide-field optical/infrared imager in space, such as the proposed 
Supernova/ Acceleration Probe (SNAP; Nugent et al. 2000) satellite. 

Fortunately, an interesting measurement of w can be made at present. 
The current values of flu from many methods (most recently WMAP: 0.27; 
Spergel et al. 2003) make an excellent substitute for those expensive SNe 
at z > 1.0. Figure 16 shows that a SN la sample with a maximum rcdshift 
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of z = 0.8, combined with the current 10% error on £2mi wm do as well 
as a SN la sample at much higher redshifts. Within a few years, the Sloan 
Digital Sky Survey and WMAP will solidify the estimate of Qm and sharpen 
w further. 
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Figure 16 (left): Constraints on fl x and w from SN data sets collected at 
z = 0.2 (solid lines), z — 0.7 (dashed lines), and z — 1.6 (dash-dot lines). 
The shaded area indicates how an independent estimate of flu with a 10% 
error can help constrain w. 

Figure 17 (right): Expected constraints on w with the desired final ESSENCE 
data set of 200 SNe la, 30 of which (in the redshift range 0.6 < z < 0.8) are to 
be observed with HST. The thin lines are for SNe alone while the thick lines 
assume an uncertainty in Qm of 7%. The final ESSENCE data will constrain 
the value of w to ~10%. 

Both the SCP and the HZT are involved in multi-year programs to dis- 
cover and monitor hundreds of SNe la for the purpose of measuring w. For 
example, the HZT's project, ESSENCE (Equation of State: SupErNovae trace 
Cosmic Expansion), is designed to discover 200 SNe la evenly distributed in 
the 0.2 < z < 0.8 range. The CTIO 4-m telescope and mosaic camera are 
being used to find and follow the SNe by imaging on every other dark night 
for several consecutive months of the year. Keck and other large telescopes 
are being used to get the SN spectra and redshifts. Project ESSENCE will 
eventually provide an estimate of w to an accuracy of ~10% (Figure 17). 

Farther in the future, large numbers of SNe la to be found by the SNAP 
satellite and the Large-area Synoptic Survey Telescope (the "Dark Matter 
Telescope" ; Tyson & Angel 2001) could reveal whether the value of w depends 
on redshift, and hence should give additional constraints on the nature of the 
dark energy. High-redshift surveys of galaxies such as DEEP2 (Davis et al. 
2001), as well as space-based missions to map the CMB (Planck), should pro- 
vide additional evidence for (or against) A. Observational cosmology promises 
to remain exciting for quite some time! 
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